Exotic Gapless Mott Insulators of Bosons on Multi-Leg Ladders 
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We present evidence for an exotic gapless insulating phase of hard-core bosons on multi-leg ladders 
with a density commensurate with the number of legs. In particular, we study in detail a model of 
bosons moving with direct hopping and frustrating ring exchange on a 3- leg ladder at v — 1/3 filling. 
For sufficiently large ring exchange, the system is insulating along the ladder but has two gapless 
modes and power law transverse density correlations at incommensurate wave vectors. We propose 
a determinantal wave function for this phase and find excellent comparison between variational 
Monte Carlo and density matrix renormalization group calculations on the model Hamiltonian, 
thus providing strong evidence for the existence of this exotic phase. Finally, we discuss extensions 
of our results to other iV-leg systems and to N- layer two-dimensional structures. 

PACS numbers: 71.10.Hf, 71.10.Pm, 75.10.Jm 



Creating a complete catalogue of the possible zero- 
temperature phases of interacting bosons in low- 
dimensional lattice systems is a fundamental open prob- 
lem in condensed matter physics [1 j. The last few decades 
have seen a number of proposals for new entries in this 
catalogue: the gamut of phases now runs from the well- 
studied "solids" and "superfluids" to exotic possibilities 
such as variants of a "Bose metal" [2H5]. Equally chal- 
lenging is the identification of parent Hamiltonians that 
realize each of the entries as a ground state. This un- 
derstanding will clearly aid in the search for novel many- 
body physics in, for example, synthesized materials [6] 
and ultracold atomic gases in optical lattices [7]. 

Weakly interacting bosons generally condense forming 
a superfluid at low temperatures. However, for strong in- 
teractions at noninteger commensurate densities, bosons 
naturally break translational symmetries of the lattice, 
crystallizing into a solid. It is also possible for these 
two orders to coexist, resulting in the formation of a 
"supersolid" for the bosons. At integer densities, there 
is also the possibility of a featureless "Mott insulator" 
that has a nondegenerate ground state with a gap to 
all excitations. These phases have been demonstrated 
to exist convincingly in numerical simulations of lattice 
models [8UTU] . States that break time reversal symmetry 
("chiral spin liquid" [11]) or lattice rotational symmetry 
("nematic" [12]) have also been proposed. A possibility 
distinct from all of these symmetry-breaking phases is a 
Bose metal, i.e., a critical quantum liquid of bosons that 
does not break any symmetries yet has many gapless ex- 
citations. Such Bose metal phases have been studied in 
detail previously, wherein the systems studied were com- 
pressible with respect to the addition of bosons [3H5] . 

In this Letter, we will propose and study a unique 
phase of matter that has gapless collective excitations 



even though the phase is an incompressible insulator 
(Fig. [I]). In particular, we will construct a trial wave 
function for this phase, suitable for use with variational 
Monte Carlo (VMC) calculations, and provide detailed 
comparisons to numerical results obtained using the den- 
sity matrix renormalization group (DMRG) method to 
argue for the existence of such an exotic phase in a spe- 
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FIG. 1: (color online) Top: Illustration of the 3- leg ladder 
geometry with periodic boundary conditions in both the x and 
y directions. The hopping strength is J in the x direction and 
J± in the y direction; the ring term K moves a pair of hard- 
core bosons on diagonal corners of a plaquette to the other 
diagonal corners (if empty). Bottom: The di and 0I2 bands 
for the 3-leg system at v — 1/3 are shown for the gapless 
Mott insulator studied here. There are exactly enough bosons 
so that the 0I2 fermions fill up one band completely. The 
bold curves indicate occupied momentum states, and note 
that there are two degenerate bands at k y — ±2tv/3. 



2 



cific microscopic lattice model of bosons. We will focus 
our study here on the square lattice 3-leg ladder (see 
Fig. [T]) at boson filling factor ^ = 1/3 as a prototypical 
example of the physics of interest. Even though our ulti- 
mate goal is to access novel two-dimensional (2D) phases, 
the more tractable ladder systems are very intriguing in 
their own right, and in our concluding remarks, we will 
propose clear extensions of our ideas that are aimed to- 
wards two dimensions. 

Wave function. — A wave function of bosons can be 
constructed by taking the product of two fermionic wave 
functions and evaluating them at the same coordinates 
(Gutzwiller projection): 



^&(ri, ...,r M ) = *di(ri, ...,r M ) ^d 2 ( r i> •> r 



M) 



(i) 



which corresponds in terms of operators to a "parton" 
construction = d\d\ with the constraint d\d\ = d\d2 = 
&tb@JG3. In the case of interest here, where time reversal 
symmetry is present on a quasi-one-dimensional ladder, it 
is natural to put each fermion in a Slater determinant by 
occupying one-dimensional (ID) bands. We stress that 
the d\ and d^ fermions need not have identical band filling 
configurations. 

For the 3-leg ladder system (Fig. [I]), there will be three 
bands of d\ and d^ fermions. At exactly v = 1/3 filling, 
an interesting possibility arises that is special to this com- 
mensurate density: there are enough fermions to exactly 
fill one c?2 band completely, but still have up to three 
partially filled d\ bands. The resulting projected state of 
bosons is an exotic "gapless Mott insulator" (GMI): the 
di state is a band insulator while the d\ state is a band 
metal. It is this phase that we will explore in detail in 
this Letter. 

Model. — The microscopic model for which we find the 
GMI phase consists of hard-core bosons hopping on a 
square lattice with a four-site ring-exchange interaction: 

H JK = - J^(&J&p+a + h - c - J ±- J2( b r b *+y + hx ') 

r r 

+ K^blbr+d+x+ybr+y + h.C.), (2) 
r 

where b r destroys a boson at site r. The J and J± terms 
are the usual hopping terms for bosons in the directions 
parallel and perpendicular to the chains, respectively. 
The ring term K is the interaction that makes the phase 
diagram interesting. When K > 0, which is the case 
studied here, the model has a sign problem thereby ren- 
dering quantum Monte Carlo calculations inapplicable. 

The K term causes a pair of bosons to hop without 
center-of-mass motion, as illustrated in Fig. [I] In this 
sense, it frustrates the attempt of the simple hopping 
term to cause the bosons to maximize their kinetic en- 
ergy by condensing and yet at the same time disfavors 
the formation of a simple solid. A more compelling (al- 
beit more abstract) motivation to study this model comes 
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FIG. 2: (color online) Summary of our numerical results: 
phase diagram in the K/J — J±/J plane of the 3-leg lad- 
der at v = 1/3 filling as determined by DMRG and VMC 
calculations. 



from a strong-coupling analysis of the U(l) gauge theory 
for the d\ and d<i fermions: within this framework, one 
can reverse engineer Eq. {2) as a potential parent Hamil- 
tonian for ground states of the form given by Eq. ([!]) 
(see Ref. [4 ). Indeed, it was found that taking d\ (cfe) 
to hop more easily along x (y) and then increasing this 
anisotropy effectively increases K/J. At v = 1/3 on the 
3-leg ladder, one hence expects that as the K term in- 
creases, the three d\ bands should become almost degen- 
erate, while the d<i bands should become almost disper- 
sionless and split apart. This in turn naturally causes the 
d2 fermions to go into a band insulating state and the d\ 
fermions to go into a metallic state (see Fig. [I]). This 
is exactly the situation discussed above in our descrip- 
tion of the GMI phase, and is our primary motivation to 
search the phase diagram of the model, Eq. ([2|, for the 
GMI phase. 

3-leg ladder. — We henceforth consider Eq. Q on the 3- 
leg ladder at filling factor v = 1/3. When J ^> if, J_l, the 
ground state is a superfluid (SF) with quasi-long-range 
order. On the other hand, when J± ^> K, J, one boson 
on each rung goes into a zero ^-momentum state, and 
the system is a conventional rung Mott insulator with a 
unique ground state and a gap to all excitations. The 
region where we expect the GMI phase is when the K 
term dominates. 

As we demonstrate below, we do in fact have exten- 
sive numerical evidence for the existence of an exotic 
GMI phase in the expected region of the phase diagram. 
Specifically, in addition to the SF and fully gapped rung 
Mott state, we find a gapless Mott insulator with two 
ID gapless modes. Our results are summarized in Fig. [2] 
We first searched the phase diagram using the variational 
wave functions: projected d\ and tfe bands with different 
band fillings (along with variational exponents on the two 
Slater determinants to give additional variational free- 
dom [5]) and Jastrow wave functions for the superfluid. 
We then also scanned the phase diagram with DMRG 
at the points marked as circles, triangles, and squares in 
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Fig. |2j keeping between 2000 and 4000 states per block 
to ensure converged results. Figure [2] shows the phase 
diagram of K/J vs. Jj_/J with colored regions (roughly) 
distinguishing the various phases. These regions gener- 
ally coincide with the VMC state that minimizes the en- 
ergy at each point, although some adjustments have been 
made to incorporate strong conclusions from the DMRG. 
The GMI phase has three partially occupied d\ bands and 
one fully filled d^ band — this is the phase of main inter- 
est. The rung Mott phase is modeled well in the VMC 
calculations by a projected wave function in which both 
the d\ and ^2 fermions fully fill the k y = band. There is 
also a SF phase and two unidentified regions (see below) 
marked as white. 

We shall focus here on establishing the existence and 
properties of the GMI phase; these points are marked 
by squares (green) in Fig. [5J A comparison between 
VMC and DMRG calculations of the Fourier transformed 
single-boson Green's function and density-density corre- 
lation function, for a representative point in the GMI 
phase, is shown in Fig.[3]for a 3-leg system of length L x = 
24. Specifically, we plot the boson momentum distribu- 
tion n 6 (q) = l/(^L x )^ r y( b t b r')e ici ' ir ~ r,) and density- 
density structure factor Db(q) = 1/(3L X ) ^ r r , ([n r — 

v][n r > — ^])e 2q '( r_r \ where n r = 6j6 r . 

In the proposed GMI phase, since the cfe fermion 
is gapped, there is a gap to any excitation that car- 
ries nonzero boson number. This implies a single-boson 
Green's function that decays exponentially in real space 
and thus a boson momentum distribution with no singu- 
lar features; the latter is clearly demonstrated in the top 
panel of Fig. [3] in both the DMRG and VMC. Because of 
the fully filled di band, the VMC state contains exactly 
one boson on each rung and hence gives completely flat 
n b(Qx,Qy = 0, ±27r/3), while the true ground state (as 
obtained by DMRG) does contain some short-distance 
fluctuations. Within the GMI, the q = mode occupa- 
tion, n&(q = 0), does not increase with increasing system 
size in the DMRG, and a finite-size study of the one- and 
two-boson gaps indicates that these gaps attain rather 
large, finite values in the L x — >> oo limit [13]. All of these 
findings are consistent with the insulating properties of 
the GMI phase. 

Measurement of the density-density correlation func- 
tion gives compelling evidence for the presence of the 
GMI phase, in particular, its gapless nature and incom- 
pressibility. In the VMC wave function, (q) contains 
contributions from gapless particle-hole excitations aris- 
ing from partially filled bands of either fermion species 
(di only in the case of the GMI) and hence has singular 
structure at various momenta [4]. Indeed, the DMRG 
calculation of the density-density structure factor dis- 
plays singular features that are impressively well repro- 
duced by the corresponding VMC calculation in which 
the lowest energy state has a band filling configuration 
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FIG. 3: (color online) Comparison of DMRG and VMC cal- 
culations of the boson momentum distribution rib(q) (top 
panel) and density-density structure factor Db(q) (middle 
panel; inset is a zoom of q y = curves) in the GMI phase 
of the 3-leg ladder at the point K/J — 2.7 and J±/ J = 1.0. 
The lack of singular features in rib(q) and q% dependence of 
Db(qx,q y = 0) around q x = 0, both highlight the insulat- 
ing nature of the GMI phase. However, the density-density 
structure factor still "sees" the gaplessness of the phase due 
to the partially filled d\ bands |4], hence the singular features 
in Db(q x ,q y = ±27r/3). These singularities can be identi- 
fied with "2/cf" wave vectors in the d\ fermion's band filling 
configuration (bottom panel) in which the k y = band con- 
tains 12 fermions and the k y — ±27r/3 bands each contain 6 
fermions. Of course, the di fermion configuration (not shown) 
consists only of one fully filled k y = band. 



shown in the bottom panel of Fig. [3] In this case, the 
main singularities in ^(q) at (q x ,q y ) — (±37r/4, ±27r/3) 
can be identified with the "2k p" wave vectors associ- 
ated with moving a d\ fermion from the k y = band 
to the k y = ±27r/3 bands (going across the bands); the 
other singularities can similarly be identified with other 
"2kp" moves that are allowed. In particular, singular- 
ities at q y = are suppressed by d^ the VMC wave 
function has exactly one boson on each rung which pro- 
duces Db(q x ,q y = 0) = 0. The DMRG calculation, on 
the other hand, gives small but finite Db(q x ,q y = 0) with 
dependence around q x = 0. This result directly im- 
plies that the system is an incompressible insulator and 
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FIG. 4: (color online) Effective central charge c as inferred 
from scaling of the entanglement entropy versus subsystem 
size S(x). For the 3-leg ladder, we fit the von Neumann en- 
tanglement entropy of leftmost contiguous blocks of size 3x 
to the scaling form S(x) = (c/3) In [(L x /tt) sin (ttx/L x )] + A 
to extract c [T5J- In the main plot, c is plotted versus K/J at 
J±/J = 1.0 for a system of length L x = 24; the jump from 
c ~ 1 to c ~ 2 at the SF-GMI boundary is striking. In the 
inset, we plot S(x) in both the SF (K/J = 1.0) and GMI 
(K/J — 2.7) phases; the points are the DMRG data, while 
the curves are the fits to the scaling form. 



rules out the possibility of a superfluid, either conven- 
tional or paired. These results are general for the entire 
GMI phase in Fig. |2j except that the locations of the 
singularities in Db(q x ,q y — ±27r/3) can be changed by 
tuning K/J and/or J±/J to alter the resulting d\ band 
filling configuration. 

A crucial aspect of the GMI phase is the expectation 
that it has two ID gapless modes at low energy: at the 
mean-field level, Fig. [I] there are naively three ID gap- 
less modes corresponding to the three partially filled d\ 
bands. We must, however, enforce a constraint such that 
d\(r)di(r) = dl(r)d2(r) to realize the physical bosons. 
This is accomplished by means of a gauge field along 
with corresponding (opposite) gauge charges attached to 
each flavor of parton. In the strong-coupling limit, one 
linear combination of the original gapless modes is ren- 
dered massive leaving only two gapless modes [5] [14] . We 
have calculated the entanglement entropy of the ground 
state wave function with DMRG and extracted from it 
the effective central charge p~5j [16] , as shown in Fig. [4] 
In the superfluid phase, we find that the central charge 
c ~ 1 as expected for a single gapless boson density mode 
in this phase. As one approaches the transition into the 
GMI phase, one finds a sharp increase to c ~ 2, providing 
strong evidence that the GMI phase has two ID gapless 
modes. That is, as elucidated by Fig. [4] even though the 
GMI is an insulator, it actually contains more entangle- 
ment and more ID gapless modes than the SF. 

In addition to the phases mentioned above (SF, rung 
Mott, and exotic GMI), there are two unidentified 
(white) regions in the phase diagram of Fig. [2j a small 
region at intermediate K/J and very small J±/J and a 



larger region above the GMI. The corresponding DMRG 
points are marked as question marks and squares (green) 
with a question mark, respectively. In the latter region, 
VMC favors a GMI ground state, but the features of the 
GMI in the DMRG fail to persist, e.g., the singularities in 
Db(q x ,Qy — ±27r/3) smoothen and the calculated central 
charge drops below c c± 2. An analysis of the spectral gap 
versus L x using exact diagonalization and DMRG reveals 
that the system is either gapless or has a very small gap 
in this region [13], without any obvious ordering. 

Extensions. — A natural extension of our findings is to 
the general TV-leg ladder at v = m/N, with integer m. 
In such systems, a GMI phase analogous to the one dis- 
cussed above is a clear ground state candidate. We leave 
for future work a numerical investigation using DMRG 
to assess the stability of GMI phases in these systems; al- 
though we can report that on the 2-leg ladder at v = 1/2 
the GMI is likely unstable in our model [14|, a result 
which partially motivated the present study. 

Another interesting extension of our results is to struc- 
tures that consist of N layers of 2D systems at a general 
filling of v = m/N e.g., a three-layer system at v = 1/3. 
In such a system, the d<i fermions can exactly fill m of the 
TV 2D bands resulting in a gapless Mott insulating state. 
Theoretically, this state is stable and can exist in two 
dimensions. We expect that an appropriate model that 
realizes this physics would have substantial ring exchange 
on the plaquettes between the layers, as this interaction 
will give the necessary disparity in hopping to keep the d\ 
system metallic while forcing the d<i system into a band 
insulating state. Thus, we are hopeful that an appropri- 
ate JK model will realize this quasi-2D analogue. 
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